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Abstract 

1 calculate the dispersion relation for quarks of mass m and momentum A; in a quark 
gluon plasma at temperature T, in the limit + k"^ ^ {qTY-, where g is the strong 
coupling constant. 1 find three contributions to the dispersion relation: one that 
depends on T but not m or /c, one that depends on m and T but not fc, and third 
contribution that depends on all three (and is opposite in sign to the other two). 
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Recently, there has been much interest in calculating the production rates for both 
massless and massive quarks in ultrarelativistic nuclear collisions Unfortunately, 
there have been no rigorous thermal field theory calculations to compare to 
approximate results, probably because the "naive" production rates diverge for 
massless quarks. The production rates can be regularized by summing over hard 
thermal loops 0, ^, following the treatment of Braaten and Pisarski [^. However, this 
is a difficult calculation, and has not been done so far. 

The closest approach to a full calculation has been to treat the thermal quark masses 
for massless quarks as if they were bare masses, in order to regulate the divergent 
production rates However, the thermal corrections for massive quarks are not 
included, and these can be important when m ~ T. There have been several papers in 
which the thermal dispersion relation was calculated for massive quarks |lT0| - p!5|] , but 



these have concentrated on the low-frequency dispersion relation. Also, no calculations 
have been done for systems out of chemical equilibrium, while in the early stages of an 
ultrarelativistic nuclear collision the quarks (and anti-quarks) have densities below 
their chemical equilibrium values. 

Because the main interest here is calculating effective masses to regulate production 
rates, I calculate the dispersion relation in the high-frequency limit. This is because the 
divergences in the rates occur at small values of t or u, but at arbitrary values of s, so 
the quarks involved have momenta of order (T^ + niTY^'^ rather than of order gT. I 
calculate the thermal masses for arbitrary bare mass, following the calculation of 
Petitgirard |]12| and working to order g^, to provide a smooth transition from the 



production of massless quarks to that of massive quarks. 

I begin with a quark gluon plasma of temperature T with strong coupling constant 
g. To one- loop order, the quark self-energy ^{K) can be written as 

= ^-cm, (1) 

where 

a = ^ [Tr (^ReS) -cjTr(^ReS)] , (2) 

b = ^ \{uj^ - A;2)Tr (^ReS) - loTt (|fReS)l , (3) 

c = -^tr(ReS). (4) 
4m ^ ^ 

Here the quark four-momentum K = {u, k), u is the matter four- velocity and m is the 
bare quark mass. The Lorentz-invariant denominator of the quark propagator is then 

Tr - m - S)^ = (1 + a)^K^ + 2(1 + a)bK^u^ + 6^ - (1 - c)^m^ (5) 

and its poles give the quark dispersion relation. 



The dispersion relation in the plasma rest frame is thus [12|, 



l + a)iu + b = J{l + ayk^ + {l + cym^. (6) 
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Because E = 0{g^), I rewrite the dispersion relation (P) as 



2h 

lk^ + m'^ + 2(c - a)m2 (k^ + m2) + 0{g'') 



UJ 



for UJ > k ^ gT. I then use the relation 



UJ 



to obtain 
where 



UJ 



2 = p + ^2 ^ ^^g2^ 

Vk^ + m^ + 6 + 0{g^) 
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5 = 2{c- a)m^ - 2buj = ^tr (ReS) + ^Tr (ITReS 
The self-energy projections have been calculated by Petitgirard |T2| 
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where Ci? is the quadratic Casimir invariant of the quark representation, nB{F) is the 
Bose (Fermi) distribution function, E = (p"^ + m^Y^'^, 

[K^ - m2 + 2p{uj + A;)] [K^ - - 2p{uj - k)] 
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[K"^ - m2 + 2p{uj - k)] [K^ - - 2p{uj + k)] 
[K^ + 771^ + 2{Euj + pk)] [K^ + - 2{Euj - pk)] 



(13) 



(14) 



[K^ + m2 + 2{Euj - pk)] [K^ + - 2{Euj + pk)] 

As the logarithms are already multiplied by g'^, and I am interested only in the lowest 
order corrections to the dispersion relation, I take K"^ = m? in the logarithms, which is 
equivalent to using the identity (||), obtaining the simpler forms 

Lb{p) = 0, 

p — k 
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The thermal correction to the bare mass is then 
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It is clear from inspection that g > and r < 0. 

For regulation of production rates in the early stages of ultrarelativistic nuclear 
collisions, the most interesting case is < 0. The case /z > is discussed thoroughly in 
Refs. |Tl[] and |T5|. For simplicity, I calculate the average value of r, 

d^k ,m k. 



r a k, ,m k , , , 
/ r — , - ni^H 
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where cu = {K"^ + m^)^/^. Using r for regulation of divergences is probably the most 
sensible procedure, in the absence of a full calculation with effective thermal 
propagators and vertices. 

The most useful case for present simulations is when the gluons are chemically 
equilibrated but the quark densities are very low (— /i/T ^ 1), as these are the 
conditions under which most quark production occurs in an ultra-relativistic nuclear 
collision. In this case, 
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The integrals can be evaluated in closed form in the low- and high-mass limits: 
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They can also be easily evaluated numerically to high accuracy. I was unable to find 
accurate interpolation formulas, probably because of the non-analytic behavior of goo 
and Too at m = 0. 

I thank T. Altherr for useful discussions. This material is based upon work 
supported by the North Atlantic Treaty Organization under a Grant awarded in 1991. 
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